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^ ■ Abstract 

^h \ The ground state properties including radii, density distribution and one 

neutron separation energy for C, N, O and F isotopes up to the neutron 
drip line are systematically studied by the fully self-consistent microscopic 



X 



Relativistic Continuum Hartree-Bogoliubov (RCHB) theory. With the proton 



density distribution thus obtained, the charge-changing cross sections for C, 
N, O and F isotopes are calculated using the Glauber model. Good agreement 
with the data has been achieved. The charge changing cross sections change 
only slightly with the neutron number except for proton-rich nuclei. Similar 
trends of variations of proton radii and of charge changing cross sections for 
each isotope chain is observed which implies that the proton density plays 
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important role in determining the charge-changing cross sections. 
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Recent progresses in the accelerator and detection techniques all around the world have 
made it possible to produce and study the nuclei far away from the stability line — so called 
"EXOTIC NUCLEI" . Based on the measurement of interaction cross section with radioactive 
beams at relativistic energy, novel and entirely unexpected features has appeared: e.g., the 
neutron halo and skin as the rapid increase in the measured interaction cross-sections in the 
neutron-rich light nuclei |l|||]. 

Systematic investigation of interaction cross sections for an isotope chain or an isotone 
chain can provide a good opportunity to study the density distributions over a wide range 
of isospin [§,£|]. However the contribution from proton and neutron are coupled in the 
measurement of interaction cross section. To draw conclusion on the differences in proton 
and neutron density distributions definitely, a combined analysis of the interaction cross 
section and other experiment on either proton or neutron alone are necessary. 

The charge-changing cross section which is the cross section for all processes which result 
in a change of the atomic number for the projectile can provide good opportunity for this 
purpose. In Ref. ||, the total charge-changing cross section a cc for the light stable and 
neutron-rich nuclei at relativistic energy on a carbon target were measured. We will study 
cr cc theoretically by using the fully self-consistent and microscopic relativistic continuum 
Hartree-Bogoliubov (RCHB) theory and the Glauber Model in the present letter. 

The RCHB theory |6|-[8|], which is an extension of the relativistic mean field (RMF) 
PHTTf and the Bogoliubov transformation in the coordinate representation, can describe 
satisfactorily the ground state properties for nuclei both near and far from the /^-stability 
line and from light to heavy or super heavy elements, as well as for the understanding of 
pseudo-spin symmetry in finite nuclei |T2|-[15|. A remarkable success of the RCHB theory 
is the self-consistent reproduction of the halo in n Li [0] and the prediction of giant halo 
||. In combination with the Glauber model, the RCHB theory successfully reproduces the 
interaction cross section in Na isotopes ||. These successes encourage us to apply the RCHB 
theory to calculate the charge changing cross section of the C, N, O, F isotopes (ranging 
from the /3-stability line to the neutron drip line) on the target of 12 C reported in Ref. ||. 



With the density distribution provided by RCHB theory the total charge- changing cross 
section can be calculated based on the Glauber model and compared with the data directly 
H, as was done in Ref. 0] for the interaction cross section. Since the theory used here is fully 
microscopic and basically parameter free, we hope it provide us more reliable information 
on both the proton and neutron distribution. 

The ground state properties of C, N, O and F isotopes up to neutron drip line are studied 
first, including single neutron separation energies, density distributions and radii. Then the 
total charge- changing cross sections will be calculated from the Glauber model with the 
density obtained from RCHB calculations. 

The basic ansatz of the RMF theory is a Lagrangian density whereby nucleons are de- 
scribed as Dirac particles which interact via the exchange of various mesons (the scalar sigma 
(a), vector omega (u) and iso- vector vector rho (p)) and the photon. The a and u meson 
provide the attractive and repulsive part of the nucleon-nucleon force, respectively. The nec- 
essary isospin asymmetry is provided by the p meson. The scalar sigma meson moves in a 
self-interacting field having cubic and quadratic terms with strengths gi and g% respectively. 
The Lagrangian then consists of the free baryon and meson parts and the interaction part 
with minimal coupling, together with the nucleon mass M and ra a (g a ), % (gS), and m p 
(g p ) the masses (coupling constants) of the respective mesons: 

C = $(ip - M)ip + Id^ad^a - U(a) - \Vt^ u 

+»^ - \R^R^ + \m 2 p p^ - \F pv F^ (1) 

-g^ipaip - g^ipipip - g p ippTip - eipflip 

For the proper treatment of the pairing correlations and for correct description of the 
scattering of Cooper pairs into the continuum in a self-consistent way, one needs to extend 
the present relativistic mean-field theory to the RCHB 
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Ek is the quasi-particle energy, the coefficients £4(r) and Vfe(r) are four-dimensional Dirac 
spinors, and h is the usual Dirac Hamiltonian 

h=[a-p + V(r)+(3(M + S(r))}, (3) 

with the vector and scalar potentials calculated from: 

V(r) = g^(r) + g p f>{v)r + ~e(l - r 3 )^r(r), 

2 (4) 

S(t) = g ff a{r). 

The chemical potential A is adjusted to the proper particle number. The meson fields are 
determined as usual in a self-consistent way from the Klein Gordon equations in no-sea- 
approximation. 

The pairing potential A in Eq. (|2]) is given by 

Aa. = ^EKX^ (5) 

cd 

It is obtained from the pairing tensor k = U*V T and the one- meson exchange interaction 
Vabcd i n the pp-channel. As in Ref. [§-|| V^ cd in Eq. (|5|) is the density dependent two-body 
force of zero range: 

V(r l7 r 2 ) = ^(l + n^-r^l-pM/po). (6) 

The ground state |\&) of the even particle system is defined as the vacuum with respect 
to the quasi-particle: (3 V \^) = 0, |\l/) = Tiis /3u\— ), where |— ) is the bare vacuum. For odd 
system, the ground state can be correspondingly written as: | \I/) M = ffi I\u^n/3u\— ), where 
\x is the level which is blocked. The exchange of the quasiparticle creation operator fy with 
the corresponding annihilation operator (3^ means the replacement of the column \x in the 
U and V matrices by the corresponding column in the matrices V*, U* ||16| . 

The RCHB equations (0) for zero range pairing forces are a set of four coupled differential 
equations for the quasi-particle Dirac spinors U(r) and V(r). They are solved by the shooting 
method in a self-consistent way as 0. The detailed formalism and numerical techniques 
of the RCHB theory can be found in Ref. and the references therein. In the present 



calculations, we follow the procedures in Ref. @J|,|] and solve the RCHB equations in a 
box with the size R = 20 fm and a step size of 0.1 fm. The parameter set NL-SH |T7| is 
used, which aims at describing both the stable and exotic nuclei. The density dependent 
5-force in the pairing channel with po = 0.152 fm -3 is used and its strength Vq is fixed by 
the Gogny force as in Ref. ||. The contribution from continua is restricted within a cut-off 
energy E cut ~ 120MeV. 

Systematic calculations with RCHB theory has been carried out for the C, N, O and 
F isotopes. The one neutron separation energies S n predicted by RCHB and their experi- 
mental counterparts |18J for the nuclei 11_22 C, 13_24 N, 15 ~ 26 and 17_25 F are given in Fig.|l| 
as open and solid circles respectively. For carbon isotopes, the theoretical one neutron sep- 
aration energies for n - 18 > 20 > 22 c are j n agreement with the data. The calculated S n is less 
than (-0.003 MeV) for the odd- A nucleus 19 C which is bound from experiment. While for 
the experimentally unbound nucleus 21 C, the predicted value of S n is positive. From the 
neutron-deficient side to the neutron drip line, excellent agreement has been achieved for the 
nitrogen isotopes. Just as the other relativistic mean field approaches, the RCHB calcula- 
tions overestimate binding for 25 > 26 which are unstable in experiment. For fluorine isotopes, 
the S n in 17 ' 26 ~ 29 F are overestimated in contrast with the underestimated one in 18 F. The 
neutron drip line nucleus is predicted as 30 F. In general, the RCHB theory reproduces the 
S n data well considering that this is a microscopic and almost parameter free model. There 
are some discrepancies between the calculations and the empirical values for some of the 
studied isotopes. This may be due to deformation effects, which has been neglected here. 

The proton density distributions predicted by RCHB for the nuclei 1 o- 22 c, 12 ~ 24 N, 14_26 
and 16_25 F are given in Fig.|2| in logarithm scale. The change in the density distributions 
for each isotopes chain in Fig.|2] occurs only at the tail or in the center part as the proton 
number is constant. Because the density must be multiplied by a factor 4iir 2 before the 
integration in order to give proton number or radii, the large change in the center does not 
matter very much. What important is the density distribution in the tail part. Compared 
with the neutron-rich isotopes, the proton distribution with less N has higher density at the 



center, lower density in the middle part ( 2.5 < r < 4.5fm ), a larger tail in the outer part ( 
r > 4.5 fm ) which gives rise to the increase of r p and a cc for the proton rich nuclei as will 
be seen in the following. 

The neutron and proton rms radii predicted by RCHB for the nuclei 10_22 C 5 12 ~ 24 N, 14_26 
and 16_25 F are given in Fig.^. The neutron radii for nuclei in each isotope chain increase 
steadily. While the corresponding proton radii remains almost constant with neutron number 
for nuclei in each isotope chain except for the proton rich ones. 

To compare the charge-changing cross sections a cc directly with experimental measured 
values, the densities p p (r) of the target 12 C and the C, N, O and F isotopes obtained from 
RCHB ( see Fig.[| ) were used. The cross sections were calculated in Glauber model by using 
the free nucleon-nucleon cross section [|1| for the proton and neutron respectively. The total 



charge- changing cross sections a cc of the nuclei 10 - 22 C ; 12 - 24 N ; 14 - 26 o and 16_25 F on a carbon 
target at relativistic energy are given in Fig. [|. The open circles are the result of RCHB 
combined with the Glauber Model and the available experimental data || are given by solid 
circles with their error-bars. The agreement between the calculated results and measured 
ones are fine. 

The charge changing cross sections change only slightly with the neutron number except 
for proton-rich nuclei, which means the proton density plays important role in determining 
the charge-changing cross sections a cc . A gradual increase of the cross section has been 
observed towards the neutron drip line. However, the big error bars of the data can not help 
to conclude anything here yet. It is shown clearly that the RCHB theory, when combined 
with the Glauber model, can provide reliable description for not only interaction cross section 
but also charge changing cross section. From comparison of this figure and Fig|| we can 
find similar trends of variations of proton radii and of charge changing cross sections for 
each isotope chain which implies again that the proton density plays important role in 
determining the charge- changing cross sections. 

Summarizing our investigations, the ground state properties for C, N, O and F isotopes 
have been systematically studied with a microscopic model — the RCHB theory, where the 



pairing and blocking effect have been treated self-consistently. The calculated one neutron 
separation energies S n are in good agreement with the experimental values available with 
some exceptions due to deformation effects which is not included in the present study. A 
Glauber model calculation for the total charge- changing cross section has been carried out 
with the density obtained from the RCHB theory. A good agreement was obtained with 
the measured cross sections for 12 C as a target. Another important conclusion here is that, 
contrary to the usual impression, the proton density distribution is less sensitive to the 
proton and neutron ratio. Instead it is almost unchanged from stability to the neutron 
drip-line. The influence of the deformation, which is neglected in the present investigation, 
is also interesting to us, more extensive study by extending the present study to deformed 
cases are in progress. 

This work was partly supported by the Major State Basic Research Development Pro- 
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of China under Grant No. 10025522, 19847002 and 19935030. 
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FIG. 1. The one neutron separation energies S Q for the nuclei n 22 C, 13 24 N, 15 26 and 17 25 F 
by RCHB theory (open circles) and their experimental counterparts (solid circles). 
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FIG. 2. The proton density distributions predicted by RCHB for the nuclei 10 22 C, 12 24 N, 
14-26q an( j i6-25p - m logarithm scale. 
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FIG. 3. The neutron and proton rms radii predicted by RCHB for the nuclei 10 22 C, 

12-24 N M-26Q and 16-25 R 
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FIG. 4. The total charge-changing cross sections a cc of the nuclei 10 -22(~i 12-24^14~26q an( j 
i6-25p on a carbon target at relativistic energy. The open circles are the result of RCHB and the 
available experimental data are given by solid circles with their error-bars. 
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